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ABSTRACT 

We explore how the expulsion of gas from star-cluster forming cloud-cores due to supernova explo- 
sions affects the shape of the initial cluster mass function, that is, the mass function of star clusters 
when effects of gas expulsion are over. We demonstrate that if the radii of cluster-forming gas cores 
are roughly constant over the core mass range, as supported by observations, then more massive cores 
undergo slower gas expulsion. Therefore, for a given star formation efficiency, more massive cores 
retain a larger fraction of stars after gas expulsion. The initial cluster mass function may thus dif- 
fer from the core mass function substantially, with the final shape depending on the star formation 
efficiency. A mass-independent star formation efficiency of about 20 per cent turns a power-law core 
mass function into a bell-shaped initial cluster mass function, while mass-independent efficiencies of 
order 40 per cent preserve the shape of the core mass function. 
Subject headings: galaxies: star clusters — stars: formation — stellar dynamics 
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1. INTRODUCTION 

One of the greatest discoveries made by the Hubble 
Space Telescope is that the formation of star clusters 
with mass and compactness comparable to those of old 
globular clusters is still ongoing in violent star forming 
environments, such as starbursts and galaxy mergers. 
The young cluster mass function (CMF: the number of 
objects per logarithmic cluster mass interval, dN / dlogm) 
is often reported to be a power-law of spectral index 
a = —2 (i.e. equ ivalent to dN oc m~ 2 dm) down to 
the detection limit (jFall fc Zhang II206H : iBik et al.ll200l 
I Hunter et al.l l2003f ). Similar result s are obtained for 
the cluster luminosity function (e.g.. iMiller et al.lll997t 
IWhitmore et al. I fl999) . However, old globular cluster 
systems are found to have a CMF which is Gaussian- 
like w ith a peak mass of ~ 10 5 M Q (jKavelaars fc Hanes I 
1997). These very different CMFs between young and 
old systems present the greatest dissimilarity between 
the two populations, and perhaps the greatest barrier to 
identifying young massive clusters as analogues of young 
globular clusters. 

Yet, it is worth keeping in mind that the luminos- 
ity/mass distribution of some young star cluster systems 
does not obey a power-law. Based on high resolution 
Very Large Array observations . Dohnsoiil ()2008l ) rules out 
a power-law cluster luminosity function for the young 
star clusters formed in the starburst galaxy Henize 2- 
10. In NGC 52 53, an irregu l ar dw arf galaxy in the Cen- 
taurus Group, ICresci et al.l (|2005f ) detect a turnover at 
5 x 10 4 Mq in the mass function of the young massive star 
clusters in the central region. Last but not least, the case 
of the Antennae merger of galaxies NGC 4038 / 39 remains 
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heavily debated. While IWhitmore et al. 1 (|1999h infer 
from their HST/WFPC2 imaging data a power-law clus- 
ter luminosity funct ion with a spectral index a ~ —2.1, 
lAnders et al.l (|2007h report the first statistically robust 
detection of a turnover at My ~ —8.5 mag, correspond- 
ing to a mass of ~ 16 x 10 3 M Q . The origin of that 
discrepancy likely resides in differences in the data re- 
duction and the statistical analysis aimed at rejecting 
bright-star contamination and disentangling complete- 
ness effects from intrinsic cluster luminosity function sub- 
structures. 

These findings raise the highly interesting prospect 
that the shape of the CMF at young ages does vary from 
one galaxy to another. If it were conclusively proven that 
the shape of the initial cluster mass function (ICMF) is 
not universal among galaxies but, instead, depends on 
environmental conditions and/or on the star formation 
process, consequences would be far reaching since this 
would imply that the young CMF may help us decipher 
physical conditions prevailing in external galaxies. 

In this contribution, we report the first results of a 
project investigating the evolution of the mass function 
of star forming cores into that of bound gas-free star 
clusters which have survived expulsion of their resid- 
ual stai^lmmin£_g^s^_Asimfiar study was carried out 
bv iParmentier fc Gilmorel (|2007t ) in which they showed 
that protoglobular cloud mass functions characterized 
by a mass-scale of ~ 10 6 M evolve into the observed 
globular cluster mass function. They propose that the 
origin of the universal Gaussian mass function of old 
globular clusters thus likely resides in a protoglobular 
cloud mass-scale common among galaxies and imposed 
in the protogalactic era. Present-day star forming cores, 
however, do not show any characteristic mass-scale as 
their mass function is a power-law down to the detec- 
tion limit. Previously, Kroupa & Boily (2002) showed 
that a power-law embedded-cluster mass function can 
evolve into a mass function with a turnover near 10 5 M 
if the physics of gas removal is taken into account. Tak- 
ing advanta ge of the grid of iV-bo d y mo dels recently 
compiled bv iBaumgardt fc Kroupa I (|2007[ ). we explore 
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in greater detail how a featureless core mass function 
evolves into the ICMF as a result of gas removal, which 
we model as the supernova-driven expansion of a su- 
pershell of gas. We also explore how the ICMF re- 
sponds to model parameter variations so as to under- 
stand what may make the mass function of young star 
cluster s vary from one galaxy to another. In a com panion 
paper (|Baumgardt. Kroupa k Parmentier II2008I ). we ex- 
plore this issue from a different angle, that is, by defining 
an energy criterion, and we apply that alternative model 
to the mass function of the Milky Way old globular clus- 
ters. 

2. NOTION OF INITIAL CLUSTER MASS 
FUNCTION 

Before proceeding any further, it is worth defining 
the notion of the ICMF properly. Following the onset 
of massive star activity at an age of 0.5 to 5Myr, a 
gas-embedded cluster expels its residual star forming 
gas. As a result, the potential in which the stars reside 
changes rapidly causing the cluster to violently relax in 
an attempt to reach a new equilibrium. This violent 
relaxation takes 10 - 50 Myr during which time the 
cluster loses stars ("infant weight-loss") and may be 
completely destroyed ("infant- mortal it y") (t h is process 
has b e en studied in deta i l bv iHills I (119801): iMathieu I 
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evolution of clusters is driven by internal two-body 
relaxation p rocesses and tidal interacti ons with the host 
galaxy (see ISpitzer I I1987T ) . Following IKroupa k Boily I 
(2002) we define the ICMF as the mass function of 
clusters at the end of the violent relaxation phase, that 
is, when their age is ~ 50 Myr. 

The effect of gas expulsion depends strongly on the 
virial ratio of the cluster immediately before the gas is 
removed. If the stars and gas are in virial equilibrium 
before gas expulsion the virial ratio Q* of the stars de- 
pends entirely on the star formation efficiency (SFE, or 
e) such that = l/(2e) (where we define virial equi- 
librium to be Q = 1/2) ([Goodwin I [20081) . Formally , 
£ is the effective SFE (eSF E: see IVerschuerenl (|1990h : 
iGoodwin k Bastian I (|2006D ). that is, a measure of how 
far from virial equilibrium the cluster is at the onset 
of gas expulsion, rather than the actual efficiency with 
which the gas has formed stars. We note that even quite 
small deviations from virial equilibrium can equate to 
large differences between the eSFE and the SFE. How- 
ever, in this paper we will make the assumption that the 
stars have had sufficient time to come into virial equi- 
librium with the gas potential (ie. gas expulsion occurs 
after a few crossing times) , and so the eSFE very closely 
matches the SFE. We discuss this assumption further in 
the conclusions. 

In this case, the radius of the embedded cluster is the 
same as that of the core and the three-dimensional ve- 
locity dispersion of the stars in the embedded cluster is 
dictated by the depth of the core gravitational poten- 
tial, namely: a — (Gmc/fc) 1 ' 2 , where m c and r c are the 



mass and radius of the core, respectively, and G is the 
gravitational constant. 

A star-cluster forming 'core' of mass m c (ie. the region 
of a molecular cloud that forms a cluster) with a SFE of e 
will form a cluster of mass e x m c . However, gas expulsion 
will unbind a fraction (1 — Fbound) of the stars in that 
cluster such that the mass of the cluster at ~ 50 Myr will 
be 



TOinit = -Fbound x e x m c . 



(1) 



(Note that when inbound ~ zer0 the cluster is destroyed.) 

It is important to keep in mind that, for clusters 
whose age does not exceed, say, 10-20 Myr, the mass frac- 
tion of stars that have escaped after gas removal is still 
small, as most stars do not have a high enough veloc- 
ity to have be come spatially dissoci a ted fro m the cluster 
(see fig. 1 in [Goodwin k Bastianl (|2006f ) and fig. 1 in 
IKroupa. Aarseth k Hurley I (|200lj )) . The instantaneous 
mass of these young clusters thus follows: m c \ < £ x m c 
and, in the absence of any explicit dependence of e on m c , 
the CMF at such a young age mirrors the core mass func- 
tion. Any meaningful comparison of our model ICMFs 
therefore has to be limited to populations of clusters 
which have either re-virialised or been destroyed. 

In our Monte-Carlo simulations, we assume that the 
SFE of a core is independent of its mass. The core mass 
function is assumed to follow a power-law of spectral in- 
dex a = —2 and we describe the probability distribution 
function of the SFE as a Gaussian of mean e and stan- 
dard deviation er e , which we denote P{e) = G(e,a e ). 

Clearly, if the bound fraction of stars after gas ex- 
pulsion, Fbound, depends upon the core mass, then the 
shape of the ICMF will differ from the (power-law) core 
mass function, an effect whose study was pioneered by 
IKroupa k Boilvl §002). We now take this point fur- 
ther by combining the deta iled A-body model grid of 
iBaumgardt k Kroupa I (|2007[ ), which provides the bound 
fraction Fbound as a function of the SFE e and of the ratio 
between the gas removal time-scale tgr and the proto- 
cluster crossing-time r cross . 

3. THE GAS EXPULSION TIME-SCALE 

We model the gas expulsion from an (embedded) clus- 
ter as an outwardly e xpanding supershell whose radius r s 
varies with time t as ([Castor. McCrav k Weaverlll975D : 
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Here Eq is the energy input rate from Type II supernovae 
and p g is the mass density of the residual star forming 
gas in the core. We assume the mass density profiles of 
the core and of the newly formed stars to be alike, that 
is, p g = (1 - e)p c . 

The energy input rate from massive stars expelling the 
gas is 

• A GR x £ S n , . 

Eq = , (3) 

where Aqr is the number of supernovae expelling the 
gas, -Esn is the energy of one single supernova, and tgr 
is the gas removal time-scale. At this stage, however, we 
know neither Aqr nor tgr- Since the energy input rate 
from massive stars is approximately constant with time 
over the whole supernova phase, we thus rewrite: 



Initial cluster mass function 



3 



E = 



N S n x E SJ , 



At 



SN 



(4) 



AtgN is the supernova phase duration (i.e. the time from 
the first to the last supernova) and -/VgN is the total num- 
ber of supernovae formed in the core. We caution that 
a fraction only of these contribute to gas expulsion, that 
is, 7V S N > N G n and AtsN > t G r. 

The gas removal time-scale corresponds to the instant 
when the supershell radius is equal to the core radius 
r c , namely, r s (t — tgr) = r c . The combination of that 
constraint on tqr with eqs. [5] and 2] gives: 



tgr. 
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(5) 



While gas-embedded cluster masses vary by many or- 
ders of magnitude, their radii are remarkab l y con stant, 
always of order lpc (see table 1 in iKroupal (|2005[ )1. As 
N$n oc m c , it follows from eq. [5] that, for a given e, the 
gas removal time-scale tqr is independent of the core 
mass m c . However, the impact of gas expulsion upon a 
cluster is not governed by the absolute value of tgr, but 
by its ratio with the core crossing time r cross , ie. Fbound 
is a function of tgr/t ci . oss . The dynamical crossing-time 
of the star forming core is given by r cross = (15/irGp c ) ' 5 
which, combined with eq. [5j gives: 
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(6) 

In this equation, the number of supernovae, N$n, from 
eq. [5] is expressed as a function of the embedded cluster 
mass m ec i = e x m r , assumi ng a Kroupa stellar initial 
mass function (Kroupa 2001). We also assume that the 
up per stellar mass limi t corre lates with m ec ;, as found 
bv IWeidner fc Kroupa I (|2004f h although the hypothesis 
of a fixed value for the mass of the most massive super- 
nova progenitor (say, 60 M ) hardly affects the results 
presented in the next section. E§\ and At$ are the en- 
ergy released per supernova in units of 10 51 ergs and the 
duration of the supernova phase, in Myr, respectively. In 
what follows, we adopt E51 = 1 and At^ = 30. 

Eq. [6] shows that, in the absence of a mass-radius rela- 
tion, the gas removal time-scale, expressed in units of the 
core crossing time, explicitly depends on the core mass. 
For a given e, the deeper potential well of more massive 
cores slows down gas expulsion when compared to low- 
mass cores. Massive clusters are therefore better able 
to adjust to the new gas-depleted potential, and have a 
larger bound fraction inbound- Eq. [6] is shown graphically 
in Fig. [1] as a match between bottom and top x-axes for 
e = 0.33 and r c = lpc. In the next section, we explore 
how the different ICMFs can be produced from the same 
(power-law) core mass function. 

4. FROM THE CORE MASS FUNCTION TO THE 
INITIAL CLUSTER MASS FUNCTION 

The core mass function is modelled as a power-law of 



spectral index a = —2 (i.e. diV 



~dm c ) and the dis- 



tribution function for the core (local) SFE is a Gaussian, 
as defined in Section [2] Both distributions are sampled 
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Fig. 1. — Mass fraction -Fbound of stars remaining bound to 
the cluster at the end of the violent relaxation phase induced by 
gas expulsion. -Fbound varies with the SFE e and with the ratio 
between the gas removal time-scale tqr and the core crossing- 
time Tcross- The arrow indicates the bound fraction -Fbound an d 
the core mass m c corresponding to the ICMF turnover when the 
mean SFE is e = 0.20 and the standard deviation of the SFE 
probability distribution is <r e = 0.00 (also marked by an arrow in 
the top panel of Fig. |2j 



independently, that is, the SFE is assumed to be inde- 
pendent of the core mass. The core radius is set to be 
r c ~ 1 pc and the ratio TGR/T cross is given by eq. [H 

The fraction -Fbound of stars remaining bound to the 
cluster after gas r emoval is inferred by lin e arly i nterpo- 
lating the grid of Ba umgardt fc Kroupa I (|2007ft (their 
table 1) which describes how Fbound varies with e and 

TGR/Tcross- 

According to our gas removal model, tgr is the time- 
scale over whi ch the cluster expels the ent irety of its 
gas. However, iBaumgardt fc Kroupa I (|2007l ) model gas 
removal as an exponential decrease with time of the clus- 
ter gas content, and their gas removal time-scale (tm in 
their eq. 1) corresponds to the time when the residual 
gas has a fraction e _1 = 0.37 of its initial value. Prior 
to using their results, we multiply tm by a factor of 3 
(i.e. tgr = 3tm), so that tgr. now corresponds to a 
residual gas mass fraction of e~ 3 = 0.05, i.e., the clus- 
ter is practically devoid of gas. The corresponding grid 
is shown in Fig. [TJ In this paper, we focus on clusters 
evolving in a low-density environment/ weak tidal field 
and we thus consider the lowest half-mass radius t o tida l 
radius ratio in the grid of IBaumgardt fc Kroupa I {2007), 
i.e. Th/r t = 0.01. In stronger tidal fields the shape of the 
ICMF is likely to be more affected than what we find be- 
low, since low-mass clusters have smaller tidal radii and 
are, therefore, more easily disrupted by gas expulsion. 

Our Monte-Carlo simulations, based on eq. [T] and 
shown in Fig. [2J show the evolution of the star forming 
core mass function into the ICMF. We have considered 9 
different cases, combining three mean star formation ef- 
ficiencies (e = 0.20, 0.33, 0.40) with three standard devi- 
ations (er £ = 0.00, 0.03, 0.05). Examination of the panels 
of Fig. [5] shows that depending on the assumed distri- 
bution of SFEs G(e, a e ), markedly different ICMFs are 
produced, ranging from a power-law shape (e.g. when 
e = 0.40, regardless of cr e ) to a bell-shape (e.g. when 
e = 0.20 and a £ < 0.03). 
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Fig. 2. — How the form of the ICMF produced from a —2 power- 
law cluster mass function responds to variations in the SFE Gaus- 
sian distribution G(e, cr E ), where e and a e are the mean SFE and 
the standard deviation, respectively. Solid lines depict the original 
core mass function. 



These different behaviours stem from how large the 
variations of the bound fraction inbound with the core 
mass m c are. In order to illustrate this, let us focus on 
the top panel of Fig. [2] for which a e — 0. A SFE as 
high as e = 0.4 guarantees that Fbound is only a weak 
function of the core mass m c as 0.35 < -Fbound ^ 0.95 
over the core mass range 10 3 M Q < m c < 10 8 M Q (see 
Fig. [I]). As a result, the shape of the ICMF mirrors that 
of the core mass function and is close to a power-law of 
spectral index a = — 2. 

In the case of a SFE of e = 0.33, the same core mass 
range corresponds to an increase of the bound fraction 
with the core mass of almost two orders of magnitude (i.e. 
0.01 < Abound < 0.90), making the ICMF significantly 
shallower than the core mass function. It is worth keep- 
ing in mind that the detailed form of this ICMF is uncer- 
tain, as the determination of the bound fraction inbound 
depends on the fine details of the TV-body modelling of 



gas expulsion when inbound J$ 0.1. As for e = 0.33, 
Abound ^ 0.01-0.04 when \og(T GR /T cross ) < -0.5 (see 
Fig-DJ. If the bound fraction were larger over that range 
of gas removal time-scale to crossing time ratio, the shape 
of the ICMF would be closer to that of the core mass 
function, while if Fbound were zero before increasing at, 
say, m core ~ 10 5 — 3 x 10 5 M Q , this would result in a 
bell-shaped ICMF (see the case of e = 0.20 below). 

The core mass function is most affected when e = 0.20. 
Fig. Q] shows that the formation of a bound gas-free 
star cluster requires its parent core to be more massive 
than 10 6 M Q , since the bound fraction is zero otherwise 
(note that the top x-axis is rightward-shifted by ~ 0.1 
in log(TG_R/T cross ) when e = 0.20, see eq. Owing to 
the wide variations of Ff, OU nd over the core mass range 
(0 ^ Abound ^ 0.85), the transformation of the core mass 
function into the ICMF is, in this case, heavily core mass 
dependent: the featureless power-law core mass func- 
tion evolves into a bell-shaped ICMF. The time-scale for 
this evolution is that required for the exposed cluster 
to get back into virial equilibrium following gas expul- 
sion, that is, 50Myr. Examination of th e model grid 
computed bv lBaumgardt fe Kroupal (|2007l ) for the SFE 
and the gas removal time-scale of relevance (e = 0.20 
and tgr > T cross ) shows that the end of the violent re- 
laxation phase occurs at about that age and the cluster 
mass is then as defined in eq. [T] 

As indicated by the arrows in Fig. [T] and in the top 
panel of Fig. [2j the turnover location is determined by 
the core mass at which the bound fraction is a few per 
cent. For instance, a core mass of 2 x 10 6 M Q with e = 
0.20 leads to log(r Gfl /r CI . oss ) = 0.13 and Fbound = 0.06. 
The corresponding initial cluster mass is thus m^t = 
(0.06) x (0.20) x (2 x 10 6 M Q ) ~ 2.4 x 1O 4 M , roughly 
the mass of the ICMF turnover. 

That e = 0.20 leads to a bell-shaped ICMF directly re- 
sults from the zero bound fraction of stars for core masses 
less than ~ 10 6 M o . A SFE of e = 0.25 would result in 
the same effect (see Fig. [J). This is equivalent to creating 
a truncated power-law for the cluster forming cores, even 
though the star forming core mass function is a feature- 
less power-law down to the low-mass regime. That result 
is reminiscent of the detailed i nvestigation performed by 
iParmentier fe Gilmore I <|2007f ) who found that a system 
of cluster parent clouds devoid of low-mass objects, such 
as a power-law mass function truncated at low-mass, re- 
sults in a bell-shaped ICMF. We emphasize that, if these 
clusters were observed at an age of 50 Myr, their observed 
mass would closely match the mass predicted by our 
model. Actually, all clusters of the bell-shaped ICMF, 
irrespective of their mass, stem from gaseous progenitors 
more massive than 10 6 Mq. This implies that unbound 
stars formed in the core leave the exposed cluster with 
velocities of order 60km.s _1 and do not linger as part of 
a halo around the bound part of the cluster. As a related 
consequence, note that if the very early Milky Way disc 
went through a violent star formation phase producing 
such massive clusters, then the corresponding high veloc- 
ity dispersions l ead naturally t o the creation of a thick 
disc component (jKrouna Il2002f) . 

In the middle and bottom panels of Fig. [2] we show 
the ICMF that results from the core mass function for 
SFE distributions with variance a £ = 0.03 and 0.05 re- 
spectively. The general features of the panels are the 
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Fig. 3. — How the shape and position of the ICMF produced 
by a —2 power-law core mass function with a SFE of e = 0.20, 
responds to variations in the width of the SFE distribution, and in 
the core radius 

same as for zero variance, however there are increasing 
differences, especially at the low-mass end of the ICMF. 
These differences are due to the effect of the SFE not 
being symmetric around the mean e. For e = 0.2, low- 
mass gas-embedded clusters with s = 0.2 are completely 
destroyed (inbound ~ 0) and this is the case for any clus- 
ter selected below the mean. However, a low-mass gas- 
embedded cluster can be drawn with e — 0.3 - 0.4, in 
which case it will be able to survive and retain a sig- 
nificant fraction of its mass. For this reason, the ICMF 
flattens rather than turning-over in the bottom panel of 

Fig. m 

A key assumption underlying the above analysis is the 
hypothesis of a constant radius for all star forming cores, 
regardless of their mass. In Fig. [3J we examine the effect 
of the core radius on the form of the ICMF produced from 
a core mass function with e = 0.20 (ie. a bell-shaped 
ICMF). We explore how the cluster mass at the turnover 
responds to core radius variations. A linear scaling for 
the y-axis is adopted as it enables us to pinpoint the 
turnover location more accurately. For the sake of com- 
parison, the ICMFs obtained with a e = 0.00 and 0.03 
and with r c = 1 pc (dotted curves with crosses in top and 
middle panels of Fig. [3J respectively) are shown as solid 
lines. ICMFs depicted with open symbols correspond to 
constant core radius of r c — 0.5 pc and r c = 2pc. A 
larger core radius implies smaller core densities, thereby 
increasing the disruptive effect of gas expulsion and low- 
ering the number of bound gas-free star clusters. That 
the turnover is at higher cluster mass arises from the 
leftward-shift of the mass-scaling of the top x-axis of 
Fig. Q] when the core radius is larger (see eq. 

As stated in Section [3l observations support the as- 
sumption of an almost constant core radius. This, how- 
ever, is at variance with theoretical expectations that the 

mass -radius relation of viriali sed cores obeys r c cx rn~J 2 
(e.g. Harris k. Pudritz I Il994f) . ICMFs obtained with 
r c [pc] = 1 0~ 3 (m c /l Mq) 1 / 2 are shown as lines with filled 
symbols in Fig. [3] The normalisation is such that a core 
of mass 10 6 M Q has a radius of 1 pc. Higher normal- 
isations, corresponding to less dense cores, lead to the 
disruption of practically the whole original population of 



embedded clusters. While, for a constant core radius of 
r c = 1 pc, the bell-shape of the ICMF is preserved when 
er e is increased from 0.00 to 0.03, the virial mass-radius 
relation leads to a power-law ICMF for a e = 0.03. That 
is, the sensitivity of the ICMF shape to the width of the 
distribution function for the SFE is greater than in case 
of a constant core radius. 

5. SUMMARY 

The shape of the mass function of young star clusters 
remains much debated. Observational results are sorted 
into two distinct, often opposed, categories: power-law or 
bell-shaped mass functions. We have shown in this con- 
tribution that both shapes can arise from the same phys- 
ical process, namely, the supernova-driven gas expulsion 
from newly formed star clusters on a core mass dependent 
time-scale. Specifically, more massive cores have a higher 
gas removal time-scale to dynamical crossing-time ratio 
by virtue of their deeper potential wells (assuming the 
near constancy of core radii over the core mass range). 
For a given SFE, this results in larger bound fractions 
Abound of stars at the end of the violent relaxation which 
follows gas expulsion. The shape of the ICMF is chiefly 
governed by the probability distribution of the SFE e. 
Bell-shaped ICMFs arise when most star forming cores 
have e < 0.25. In that case, however, bound gas-free 
star clusters only form out of cores more massive than 
1O 6 M . This implies that bell-shaped ICMFs are also 
associated with gas-rich environments so as to guarantee 
that the core mass function is sampled to such a high 
mass. In contrast, power-law ICMFs mirroring the core 
mass function arise when e ~ 0.4 for most star forming 
cores. This implies that the shape of the ICMF also de- 
pends on the width a E of the SFE distribution. When 
£ < 0.25, an increasing o~ e enhances the contribution of 
cores with e ~ 0.40 and turns a bell-shaped ICMF into 
a power-law. The amplitude of that evolution is greater 
when core radii are assumed to follow the mass-radius 
relation of virialised cores than when the core radius is 
assumed to be constant. As illustrated by Figs. [2] and 
[3j that the mass function of young clusters observed in 
a variety of environments appears not to be universal is 
n ot surprisin g. 

IVesperinTI (|1998h pointed out that a power-law ICMF 
of spectral index (-2) evolves to a bell-shaped mass func- 
tion the turnover of which is located at a significantly 
lower cluster mass than what is observed for the uni- 
versal globular cluster mass f unction. For such models, 
iParmentier fc Gilmore I (|2005f ) showed that a power-law 
ICMF with a = — 2 also leads to a radial gradient in 
the peak of t he evolved mass function contrary t o what 
is observed. iBaumgardt. Kroupa fc Parmentierl (2008) 
demonstrated that these two issues do not arise if the 
ICMF is depleted in low-mass clusters compared to a 
power-law of spectral index (-2) as a result of gas expul- 
sion. In that case, it evolves within 13 Gyr into a cluster 
mass function similar to that of old globular clusters both 
in the inner and outer Galactic halo (see their fig. 4). 
While a mean SFE of e > 0.40 preserves the shape of the 
core mass function, a mean SFE not exceeding e ~ 0.25 
leads to a low- mass cluster depleted ICMF (see Fig. [2]). 
The universality of the globular cluster mass function 
may thus stem from the existence of an upper limit of, 
say, 25 per cent on the mean SFE in all protogalaxies. 
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The origin of that constraint on the mean SFE in the 
protogalactic era remains to be explained, however. 

In this preliminary study, the energy input rate from 
massive stars encompasses the contribution made by su- 
pernovae only, i.e. that of stellar winds has been ne- 
glected. Consequently, the above results apply to the sole 
case of metal-poor environments. As noted in Section 2, 
we have assumed that the effect of gas expulsion depends 
on the actual star formation efficiency - i.e. we assume 
that the stars and gas (potential) are in virial equilib- 
rium at the onset of gas expulsion. This assumption is 
reasonable if (the stellar component of) the cluster has 
had time to relax. A cluster will have been able to relax 
if the onset of gas expulsion occurs after a few crossing 
times. In a low-metallicity regime such as we are con- 
sidering in this paper, gas expulsion will be driven by 
supcrnovae rather than stellar winds and so will not be- 
gin until a few Myr after the stars have formed. Such a 
time-scale is several crossing times, and so we could ex- 
pect the eSFE to closely match the SFE. We note how- 
ever, that there may well be a mass-eSFE dependence in 



higher metallicity cluster populations that may alter our 
conclusions (which would be applicable to populations 
such as those in the Antennae merger of disc galaxies). 

In a follow-up paper, we will consider the additional 
impact of stellar winds and the case of star clusters form- 
ing out of dense environments where strong tidal fields 
prevail. 
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